Exact soliton solution and inelastic two-soliton collision in spin chain driven by a 

time-dependent magnetic field 
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We investigate dynamics of exact N-soliton trains in spin chain driven by a time-dependent magnetic 
field by means of an inverse scattering transformation. The one-soliton solution indicates obviously 
the spin precession around the magnetic field and periodic shape-variation induced by the time 
varying field as well. In terms of the general soliton solutions N-soliton interaction and particularly 
various two-soliton collisions are analyzed. The inelastic collision by which we mean the soliton shape 
change before and after collision appears generally due to the time varying field. We, moreover, show 
that complete inelastic collisions can be achieved by adjusting spectrum and field parameters. This 
may lead a potential technique of shape control of soliton. 
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I. INTRODUCTION 

Over the past three decades, an enormous amount of 
literature has appeared throughout soliton physics and 
the underlying completely integrable models. The clas- 
sical Heisenberg spin chain which exhibits both coherent 
and chaotic structures depending on the nature of the 
magnetic interactions [1-4] has attracted considerable 
attentions in nonlinear science and condensed-matter 
physics. Solitons in quasi one-dimensional magnetic sys- 
tems have already been probed experimentally by neu- 
tron inelastic scattering [5,6], nuclear magnetic resonance 
[7,8], Mossbauer linewidth measurements [9], and elec- 
tron spin resonance [10]. The corresponding theoretical 
studies are based usually on the Landau-Lifshitz equation 
[11]. The isotropic spin chain has been studied in various 
aspects [12-16] and the construction of soliton solutions 
of Landau-Lifschitz equation with an easy axis has been 
also discussed [17,18]. It is demonstrated that the in- 
verse scattering transformation [14,19-21] can be used 
to solve the Landau-Lifschitz equation for an anisotropic 
spin chain. Great efforts [22,23] have been devoted to 
construct the soliton solution which are found by means 
of the Darboux transformation [24] . The continuum spin 
chain in an external magnetic field is of great interest 
and multi-soliton solutions of Landau-Lifschitz equation 
for an isotropic spin chain have been reported [25]. Us- 
ing Darboux transformation the nonlinear dynamics of 
anisotropic Heisenberg spin chain in an external mag- 
netic field is investigated and exact soliton solutions are 
obtained [26]. Recently soliton interaction has been in- 
vestigated [16]. The main goal of this paper is to study 
the new effect of soliton-soliton interaction in spin chain 
driven by time oscillating magnetic field. We obtain ex- 
act solution of N-soliton trains in terms of an inverse 
scattering transformation. It is shown that inelastic col- 
lisions generally appear due to the time- varying field and 



the complete inelastic collisions which may lead to a in- 
teresting technique of soliton filter and switch can be 
achieved in special case. 

The outline of this paper is organized as follows: In 
Sec. II the formalism obtained by an inverse scatter- 
ing transformation is explained in detail and the general 
N-soliton solution for reflectionless case is obtained. Pre- 
cession of nonlinear spin waves in the oscillating magnetic 
field is shown in Sec. III. Sec. IV is devoted to general 
two-soliton solution and soliton collisions. Finally, Sec. 
V will give our concluding remarks. 



II. EXACT SOLUTION OF N-SOLITON TRAIN 

Our starting Hamiltonian describing the spin chain in a 
time oscillating magnetic field with an arbitrary direction 
can be written as 



H = -J 



S n ■ S n > - giiB~B (t) ■ 22 S„ 



(1) 



<n,n'> 



where S n = (S%,S%,S*) with n = 1,2, ...,7V are spin 
operators, J > is the pair interaction parameter, 
g the Lande factor and fis is the Bohr magneton, 
B (t) = B cos(uit)e is the external magnetic field with 
e = (sin#, 0, cos#) denoting the unit vector of field direc- 
tion where chain axis and direction of magnetic field are 
assumed in x-z plane. The angle 9 between direction of 
magnetic field and z-axis is arbitrary. 

The equation of motion for the spin operator on the nth 
site is 4iS n = —j-[S n ,H]. At low temperature, the spin 
can be treated as a classical vector such that S n — > S(x). 
So that the equation of motion in a continuum spin chain 
under a time-dependent magnetic field can be obtained 
as a Landau-Lifschitz type 



dt 



S x ( ° 2 
\ dx 2 



(2) 



1 



with e = gfj, B B(i) / (2J), where S(x,t) = (S x (x,t), 
S y (x,t), S z (x,t)). We set the length of the spin vec- 
tor to unit for the sake of simplicity S 2 (x,t) = 1. The 
dimcnsionless time t and coordinate x in Eq. (2) are 
scaled in unit jj and d respectively, where d denotes the 
lattice constant. 

The corresponding Lax equations for the equation of 
motion (2) are written as 



—^{x,t,X) =L(\)*(x,t,\), 
§^f(x,t,X) = M(A)f(x,*,A), 



(3) 



where A is the spectral parameter, "J (x, t, A) is eigenfunc- 
tion corresponding to A, and L and M are given in the 
form 

L = — iX(S ■ a), 

M= l -{e-<r) + «2A 2 (S • a) - A(S • <^)(^S • a). (4) 

Here a is Pauli matrix. Thus Eq. (2) can be recovered 
from the compatibility condition J^L — -§^M + [L, M] = 
0. Based on the Lax equations (3), we derive the ex- 
act iV-soliton solution by employing the inverse scat- 
tering transformation. We consider the following nat- 
ural boundary condition of initial time(t = 0), S (x) = 
(S x ,S y ,S z ) — > (sin 6*, O,cos0) as |x — > oo, namely, the 
spin vector is along the field direction. We then have the 
asymptotic form of Eq. (3) at \x\ — ► oo, 



d x E(x,X)=L (X)E(x,X), 



where 



E(x,X) = Ue 



-iXxas 



L (A) = -iXU , 



and 



cosf 



suit 



Uq I , ,, 

1 smO — cos 6 



tan ; 



tan ■ 



(5) 



(6) 



• (7) 



The Jost solutions (x, A) and (x, A) of Eq. (3) are 
defined as 

&+(x, A) — ► E (x, A) as x — > oo, 

\f r _ (x, A) — > E (x, A) ass^ — oo. 

With standard procedures, one finds the following inte- 
gral representations of the Jost solutions in terms of the 
integration kernels K and N to be determined, 

/>oo 

*+(x, A) = Ue- lXxa:i + A / dyK{x,y)Ue- lXya \ 

J X 



K(x, oo) = 0, if (x, y) = as y < x. 



(8) 



(x, A) = Ue- iXxa * + X f dyN (x, y) Ve~ iXya3 , 

N(x, -oo) = 0, N(x,y) = as y < x. (9) 

where if and TV are 2x2 matrices. Substituting (x, A) 
in Eq. (8) into Eq.(3) and noting Uo z U~ x = U , we 
obtain 

S-tT=[I-iK(x,x)U ]U [I-iK(x,x)U ]~ 1 (10) 

where I is unit matrix. It is obvious that Eq. (10) gives 
rise to a relation between kernel K and spin vector S to 
be obtained. 

The scattering data for the operator L(x, A) are the set 
s = {a(X),b(X); X n ,c n , ImX > 0, n — 1,---,N}, where 
a(A)| 2 + |fo(A)| 2 = 1, and the function a(X) can be ana- 
lytically continued to the half-plane imA > 0. The dis- 
crete eigenvalues, A„, for the operator L (x, A) are zeroes 
of a(A) such that a(A„) = (for the simplicity we con- 
sider only simple zeroes). The functions a(A) and 6(A) 
are seen to be transmission and reflection coefficients of 
the operator L respectively. The parameter c„ denotes 
the asymptotic characteristics of the eigenfunctions. 

The time-dependence of the scattering data s (t) can 
be obtained from the second Lax equation (3), 

a(X,t) = a(A,0), 

b(X, t) = exp ( -4iX 2 t - r p T b(X, 0), 



c n (0). (11) 



Jlu 

X n (t) = X n (0), 

c n (t) = exp f -4«A 2 t - i j 



where c„(0), 6(A, 0) and a(X, 0) are constants deter- 
mined by initial conditions. The Gclfand-Levitan- 
Marchcnko equation establishes a relation between the 
kernel K (x,y,t) and the scattering data s (t) and has 
the form 

K(x,y,t)U^+F 1 + ^- J X- 1 r(X)F 2 dX = 0, 

(12) 

as y > x, where r (A) = b (A) /a (A) and 



x+v) 



and 



V / „=i A ™ 
+ r if (x, z,t)u(°)j2 (*) e iX ^+*Uz, 

Jx ^ ' n=l 

F 2 = £7 ^ J ^ e* A * + A if (x, z, ( 1 ) e ^ dz - (13) 

For the reflectionless case, r (A) = 0, Eq. (12) becomes 
a set of algebraic equations and after tedious calculation 
the matrix elements of the kernel if arc obtained as 



2 



Kn(x,x,t) = cos 2 -[Bi + S 2 tan-], 

9 
K 12 (x, x, t) = cos 2 - [Bi tan - - B 2 ] . 

with 

dot [7 + G"G' + 7> T (Ctan § - CG')} 

7>1 = : — — 1, 



follows 



(14) 



Bo 



det(7 + G'G") 

det[7 + G'G" -D T (C + CG" tan f )] 
det(7 + G'G") 



1. (15) 



where C(x,t), C' (x,t), D (x) are 1 x N matrices, 
G' (x, t), G" (x, t) N x N matrices, respectively. The su- 
perscript T means the transposed matrix and the overbar 
denotes complex conjugate, 

C(x,t) n = c n (t) K 1 D{x) n , 
C'(x,t) n = c n (t)D(x) n , 
D(x) n = cxp(«A„x), 

1 



G 1 (x, t) 
G" (x, t) 



i (A n — A TO ) 
1 



D(x) n C' (x,t) m , 



D(x) n C'(x,t) m . (16) 



— i (A„ — A m ) 

Substituting Eq. (14) into Eq. (10) , we obtain the 
general form of N-soliton trains, 

S x = ^ Re { - i2K 12 [1 - iK u cos 9} 

+ [1 + - K 2 12 ] sin^}, 

S y = ^ Im { - i2K 12 [1 - iK n cos 9} 

+ [l + K 2 11 -K 2 12 ] S m9}, 

S z = ^{ [1 + \K n \ 2 - \K 12 \ 2 ]cos9 

+2 Im [K u (1 + iK^sm 9)]}. (17) 

where Ki 2 is the complex conjugate of 7Ti2. 

A = 1 -i[Ku cos 9 + K 12 sin 9} | 2 



ITCii sin 9 — K 12 cos ( 



(18) 



According to exact N-soliton solutions in Eq. (17), wc, 
generally speaking, can investigate the dynamics of soli- 
ton trains and soliton interaction. The neighboring soli- 
tons may repulse or attract each other with a force de- 
pending on their phase difference. Particularly we in the 
following shall concentrate on the analyses of one-soliton 
dynamics and two-soliton collisions which may be of more 
interest. 



III. ONE-SOLITON DYNAMICS AND SPIN 
PRECESSION WITH TIME VARYING 
AMPLITUDE 

When N = 1, from Eq. (14) and Eq. (17) we ob- 
tain the general form of the exact one-soliton solution as 



S x 
S v 



Ri 



|Ai| 4 cosh 2 9i : 

#2 



Ai| 4 cosh 2 6i' 
S z = R 3 cos 9 + Ri sin 9, 



(19) 



where 



-20i 



sin 9 



i?i = [|Ai| 4 cosh 2 6i +P 2 (al -01 cos2(9)e 
+/3 2 |Ai| 2 (2cos 2 0sin 2 $i - l)sin0 
+2 fa 2 1 Ai | (2/3i sin $i sin 2 6 + ai cos $i)e" 01 cos 9 
-2 fa ((3ie~ &1 sin 9 + | Ai | sin $i cos 9)[\X 1 \ 2 cosh 9i 
+/3i(|Ai|sin$isin6»- /?ie" ei costf) cosfl], 
R 2 = 2a 1 0l\\ 1 \e- Bl sin $i + 2fa\fa | 3 cos $i cosh 9i 
— 2/?f | Ai| (sin 9 + cos 2 9) e~ Bl cos$i, 
2/?i 2 



7*3 = 1 

i?4 = 



|Ai| 2 cosh 2 ei' 

1 



|Ai| 2 cosh 2 ei 1 
+2a\fa sin($i 



2p( cos($i - ^i)sinh6i 
- 0i ) cosh 0i ] 



(20) 



with 



@i = 2fa(x-V 1 t)-x 1 , 

V x =4ai,xi = \n[(2fa)~ 1 c 1 (0)], 

$i = 2aix - 4(a 2 - fa\)t - (u>jy 1 gfi B Bsm(ujt) - fa, 
n 1 = 2a^ 1 (a 2 - fa 2 ) + fl B , 

Qb — (2ai7) 1 gfJ-sB cos uot, (21) 

fa = arg Ai, Ai = ai + ifa is eigenvalue parameter. The 
solution (19) describes a spin precession around magnetic 
field direction characterized by four real parameters: ve- 
locity Vi, frequency fii, coordinate of the center of the 
solitary wave xi and initial phase fa. The center of soli- 
tary wave moves with a velocity V\ , while the wave depth 
and width vary periodically with time. The wave shape 
is modulated periodically by frequency fli depending on 
magnetic field. Therefore, the solution (19) can not be 
written as the form of separating variables. Amplitude 
A and phase <J>i are complicated functions of J, B, u> 
and Ai. When ol\ = fa, the frequency fli depends on 
magnetic field only, and we have f2i = CIb- If ot-i = fa 
and B = 0, the solution (19) reduces to the usual soliton 
without shape changing. Therefore, we can use magnetic 
field to adjust spin precession and the wave shape as well. 

For a special case, 9 — 0, namely the magnetic field 
is along the z-axis, S z is independent of magnetic field, 
S z = 7*3, while S x and S v precess around magnetic field 
(z-axis). The precession frequency Hi is determined by 
magnetic field. As magnetic field rotates from 9 = (z- 
axis) to 9 = n/2 (x-axis), we can find the correspondence 



3 



such that S x -» -S z , S y -> S y , S z -> S x . The three 
components of spin vector satisfy "left-hand rule" . When 
6 = 7r/2, S x is independent of magnetic field, while S y 
and S z precess around magnetic field (x-axis). These 
results show that the magnetic field results in the motion 
of the center of solitary waves along the field direction 
and the spin vector rotates around the field in any case. 



IV. TWO-SOLITON COLLISION 

When N = 2, from Eq. (14) and Eq. (17) the general 
form of the exact two-soliton solution is seen to be 

S x = Re[-i2Q 2 (l -iQ 1 cos 6) + (\ + Q\- Q 2 )sin0], 
S y = lm[i2Q 2 (l-iQiCos6) -(l + Q\- Q 2 )sin6>], 
S z - (1 + |0i| 2 - |O 2 | 2 )cos0 

+2Im[Qi(l + iQ7sinfl)]. (22) 



here 



arg Xj and A,- 



- i/3j is eigenvalue parame- 



where 



Qi 



e 



W 



4(/l-/3)/6 + (/2-/ 4 )/ 5 



+ tan-[(/ 1 -/ 3 )/ 8 + (/2-/4)/7]}, 



COS 



,2 e 



Q2 = - h) h + (/2 - h) h] tan - 

-C/T-75)/ 8 -(7i-7D/7}, (23) 

with 

fl = 1 + + XlX2<7l9 2 , /2 = 1 + |g 2 | 2 +X1X29192, 

/3 = Xi l?i| 2 + Xi<7i9 2 , .A = X 2 M 2 + X29i92, 
/5 = Ci ( <?i tan ^ - | 9 i| 2 j - Xi6<Zi<7 2 , 



/e = 6 ^2 tan - - \q 2 \ j - X2£i9i<?2, 

h = -ii (qi + ki| 2 tan0 - 6Xi9i92tan^, 

fs = ~6 f <Z 2 + M 2 tan M - £iX 2 <?i9 2 tan ^, 



Xi 



2/3iAj 



-i(Ai -A 2 ) |A!|' X2 -i(A 2 -A 1 )|A 2 |' 
VF = A/ 2 - hh, q-j = e~ e ^ , & = 2/3, lA.p 1 , (24) 

and 

6^ = 2^(3;-^)-^, 
K j =4a j ,o ;j =l n [(2/3 J )- 1 c J (0)], 

&j = 2ctjx — 4(a 2 — f3j)t — (ojJ)^ 1 g/isB sin(wt) — 0j, 



2/3 2 A 2 



ils = (2a j J) 1 g/iBB cosojt, 



tcr, j = 1, 2. The solutions (22) describe a general inelas- 
tic scattering process of two solitary waves with different 
center velocities V\ and V 2 , different shape variation fre- 
quencies fix and 0,2- Before collision, they move towards 
each other, one with velocity V\ and shape variation fre- 
quency fli, the other with V 2 and 2 . The interaction 
potential between two solitons is a complicated function 
of parameters J,B,u> and Xj. When ay = f3j, two-soliton 
shape- variation frequencies £lj(j = 1,2) are determined 
by magnetic field. In the case of B = 0, the solutions (22) 
reduce to that of the usual two-soliton with two center 
velocities while without shape change where a interest- 
ing process in the absence of magnetic field is that the 
collision can result in the interchange of amplitude Aj 
and phase = 1,2) like exactly in the case of elastic 
collision of two particles. 

In order to understand the nature of two-soliton in- 
teraction, we analyze asymptotic behavior of two-soliton 
solutions (22). Asymptotically, the two-soliton waves 
(22) can be written as a combination of two one-soliton 
waves (19) with different amplitude and phase. The for- 
mation of two-soliton waves in the corresponding limits 
x — > — oo and x — > oo is similar to that of one-soliton 
waves (19). Analysis reveals that there is an amplitude 
exchange among three components S x , S v and S z of each 
soliton during collision, which can be described by a tran- 



sition matrix T, k such that A k+ 



A 



T k , where the 



(25) 



subscript 1 = 1,2 respectively represents the first and 
the second soliton, k = x,y,z denote three components 
of each soliton, the sign ± denotes the asymptotic limits 
of the corresponding amplitude, A k± , at x — > ±oo. As a 
consequence, amplitude change of the three components 
S k of the first soliton from A k ~ to A k+ is given by square 
of transition matrices \T k \ 2 along with phase shift 5$>\ 
during collision. In a similar fashion, the three compo- 
nents S 2 of the second soliton also change amplitudes 
from A 2 ~ to A k+ with a quantity |T 2 fe | 2 . The associate 
phase shift for the second soliton is £$ 2 . We also note 
a net change of separation distance between two solitons 
by SX 12 . 

For the special case \T k | = 1, which is possible only 
when A 2 = — Ai, we have the standard elastic collision. 
For all other cases, we have the quantity \T k \ /=1, which 
corresponds to relative change among three components 
of the spin vector leading to the deformation of soliton 
shape. However, the total amplitude of individual soli- 
tons Si and S2 is conserved quantity i.e., X); l^f*! 2 i s 
constant for I = 1,2. 

It is interesting to show the inelastic collision graphi- 
cally. The general inelastic head on collision is explained 
in Fig. 1 from which it is seen that the amplitudes of Si, 
S 2 are respectively suppressed and enhanced after colli- 
sion. Fig. 2 are devoted to the complete inelastic head 
on collisions. The amplitudes of Siand S 2 are respec- 
tively suppressed after collision shown in Fig. 2a and 2b. 
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The complete inelastic overtake-collision is shown in Fig. 
3 with the amplitudes of Si and S 2 suppressed, respec- 
tively. 



V. CONCLUSION 

In terms of an inverse scattering transformation the 
exact solution of N-soliton trains in a spin chain driven 
by a time oscillating magnetic field is obtained. From the 
general solution the dynamics and soliton interactions are 
analyzed. The one-soliton solution gives rise explicitly to 
the spin precession along with the soliton shape variation 
induced by the time varying field. It is also shown that 
the time varying field leads generally to the inelastic and 
particularly the complete inelastic two-soliton collisions 
which may be useful in developing a soliton-shape control 
technique. 
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Figure caption 
Fig. 1 

Inelastic head on collision between two solitons - pro- 
files of z-component S z (x; t) of spin vector in Eq. (22) in 
spin chain under a time-dependent magnetic field show- 
ing two different dramatic scenarios of the shape chang- 
ing collision, where 9 = Ai = —0.2 + «0.45, A 2 = 
0.3 + i0.65, ci (0) = -0.2, c 2 (0) = 3.5, g^ B B/J = 0.01, 
cj = 10, V\ = —0.8, V 2 = 1.2. All quantities plotted are 
dimcnsionless. The same is in Fig. 2 and 3. 

Fig. 2 

(2a) Complete inelastic head on collision expressed 
by Eq. (22) when Si suppressed, where 9 = 0, Ai = 
-0.35 + i0.4, A 2 = 0.2 - i0.6, a (0) = 0.2, c 2 (0) = -2.5, 
9Li B B/J = 0.01, u = 10, Vx = -1.4, V 2 = 0.8. 

(2b) Complete inelastic head on collision expressed 
by Eq. (22) when S 2 suppressed, where 9 = 0, Ai = 
-0.35 - iOA, A 2 = 0.2 + i0.6, a (0) = -0.2, c 2 (0) = 2.5, 
gUsB/J = 0.01, u = 10, V x = -1.4, V 2 = 0.8. 

Fig. 3 

(3a) Complete inelastic overtake-collision expressed by 
Eq. (22) when Si suppressed, where 9 = 0, Ai = —0.55 + 
i0.4, A 2 = -0.1 - i0.45, a (0) = 0.2, c 2 (0) = -2.5, 
gH B B/J = 0.01, u = 10, Vi = -2.2, V 2 = -0.4. 

(3b) Complete inelastic overtake-collision expressed by 
Eq. (22) when Si suppressed, where 9 — 0, Ai — —0.55 — 
iOA, A 2 = -0.1 + i0.45, a (0) = -0.2, c 2 (0) = 2.5, 
g^B/J = 0.01, lo = 10,Vi = -2.2, V 2 = -0.4. 
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